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ABSTRACT 

We present Submillimter Array 1.3 mm waveband continuum and molecular 
line observations of the HH 80-81 high-mass star- forming region. The dust con- 
tinuum emission reveals two dominant peaks MM1 and MM2, and line emission 
from high-density tracers suggests the presence of another core MC. Molecular 
line emission from MM1, which harbors the exciting source of the HH 80-81 radio 
jet, yields a hot molecular core at a gas temperature of 110 K. The two younger 
cores MM2 and MC both appear to power collimated CO outflows. In particular, 
the outflow arising from MM2 exhibits a jet-like morphology and a broad velocity 
range of 190 kms -1 . The outflow contains compact and fast moving molecular 
clumps, known as "molecular bullets" first discovered in low-mass class proto- 
stellar outflows. These "bullets" cannot be locally entrained or swept up from 
the ambient gas, but are more likely ejected from the close vicinity of the cen- 
tral protostar. The discovery of this remarkable outflow manifests an episodic, 
disk-mediated accretion for massive star formation. 



Subject headings: ISM: individual (HH 80-81) - - ISM: jets and outflows - 
circumstellar matter — stars: formation — stars: early-type — masers 
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Introduction 



Over the last decade, there have been a growing number of high-angular-resolution 
studies using observations from submm/mm and cm interferometers, aimed at understand- 
ing morphological and kinematical features associated with the massive star formation pro- 
cess, such as molecular outflows, possible rotating disks or envelopes, and inward accretion 
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flows fe.gJZhang et al.lll998l:IShepherd et al.lll998l : Beuther et alJl2002t ICesaroni et alJl2005 



20091 ) . More recently, studies based on observations from 



the Submillimeter Array (SMA), a combined power of high-angular-resolution and broad 
frequency coverage, have begun to shed light on the physical and chemical diversities in 
high-mass star-forming cores at different evolutionary stages (e.g., Rathborne et al. 2008; 
Beuther et al. 2009; Zhang et al. 2009). These efforts provided a wealth of information that 
greatly advanced our knowledge of massive star formation. 

In this paper we present a high-angular-resolution study of the high-mass star-forming 
region H H 80-81, which has a b olometric luminosity of 20, 000 L & at an adopted distance of 
1.7 kpc ([Rodriguez et al.lll980t ). This region is distinguished by a 5.3 pc radio jet that has 
an openin g angle of known as the best collimated radio jet driven from a high-mass 
protostar (jMarti et al.lll993l . 1 19951 . Il998l ). Recent Spitzer 8 /zm imaging reveals a bi-conical 
outflow cavity surrounding the radio jet (lOiu et al.l I2008T) . The iet powering source has 
been detected from mid-infrared to cm wavelengths (jAspin et al. 1994 ; Stecklum et al. 1997 ; 
Marti et al. 1999 ; Gomez et al. 2003 ; Qiu et al1l2008 ). yet some basic properties of the source, 
such as physical and chemical conditions of the circumstellar matter, remain unclear. Here 
we present a detailed study of the dust continuum and molecular line emission of the radio 
jet powering source, and the discovery of an exceptional CO outflow containing fast moving 
molecular bullets. 



2. Observations and Data Reduction 

The SMA0 observations toward the phase center (R.A., decl.) J20 oo = (18 h 19 m 12f09, -20° 
47 30. 9) were undertaken with eight antennas on 2008 June 30 in the compact configura- 
tion and on 2008 September 9 in the sub-compact configuration, both under good weather 
conditions. Particularly the sub-compact data were taken with t 2 25gh z ~0.02-0.03. The 4 
GHz correlator was configured to cover rest frequencies of 219.2 to 221.2 GHz in the lower 



lr The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 
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sideband and 229.2 to 231.2 GHz in the upper sideband. In the compact configuration, we 
observed Uranus, 3c454.3, and 3c273 for bandpass calibration, Callisto for flux calibration, 
and J1911-201 to monitor the time dependent gain. In the sub-compact configuration, obser- 
vations of 3c454.3, Uranus, and J1733-130 and J1911-201 were used for the bandpass, flux, 
and gain calibrations, respectively. The data were calibrated using the IDL MIR package 
and then exported to MIRIAD for further processing. Self-calibration was performed on 
continuum and the solutions were applied to lines as well. The two data sets were jointly im- 
aged with the ROBUST parameter set to 0.5, resulting in a synthesized beam of ~3."5x2."8 
at FWHM. 



Results 



3.1. Dense Molecular Cores 



Figure [T] shows the 1.3 mm continuum map and the integrated intensity maps of lines 
detected from a variety of molecular species. The mm continuum reveals two dominant 
peaks MM1 and MM2; MM1 coincides well with the radio source that powers the HH 80-81 
radio jet and MM2 is spatially assoc iated with H 2 masers and a very weak radio source 
(IGomez et al.lll995t iMarti et al.lll999l ). In addition, emission lines from high-density tracers, 
such as CH 3 OH, CH 3 CN, H 2 13 CO, OCS, and HNCO, clearly reveals another dense core MC 
(Figures Hb-QJ), although the mm continuum map only shows a weak enhancement toward 
this core. In Figure [lb the CH3OH (8-i£-7oj)E emission toward MC is eminently bright 
with a peak flux of 6.4 Jybeam -1 . This n arrow spectral com ponent at a LSR velocity of 
13.2 kms -1 appears to be a maser emission. iKurtz et al.l (120041 ) detected a 44 GHz methanol 
maser in the vicinity o f MC, which is con sistent with our 229.8 GHz detection as both masers 
are class I transitions (jSlysh et al.l 120021 1 . Apart from the maser emission, most thermal line 
emission in Figure [T] appears brighter toward MM1. In particular lines of CH 3 OH (15^11- 
16 3) i 3 )i? and HNCO (10i )9 -9i i8 ) and those from S0 2 isotopologues are only seen in MM1. 
We will discuss these diversities in Section 14.11 

MM1 has an integrated flux of 1.2 Jy at 1.3 mm continuum, in agreement with the 
BIMA measurement of 0.97 Jy at 1.4 mm and is completely dominated by the dust emission 
(IGomez et al.l 120031 ) . The mm continuum emission from MC does not appear to be very 
centrally peaked and is blended with the MM2 emission. Therefore we cannot obtain a 
reliable estimate of the continuum flux for MC alone. The measured flux toward MM2 and 
MC amounts to 0.8 Jy and is apparently dominated by MM2. From the continuum fluxes and 
assuming thermal equilibrium between the gas and dust hence adopting the temperatures 
derived from the CH3CN emission (Section 14. ip . we obtain a gas mass of 6-34 M & for MM1 
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and a total mass of 12-65 M Q for MM2 and MC, fo r a dust opacity K225GHz~1-9-0.35 cm 2 



g 1 or an emissivity index /3~l-2 (IHildebrand 1119831 ) 



3.2. Molecular Outflows 

Figure [2] shows the high velocity CO (2-1) emission integrated over every 20.4 kms" 1 
from —85.2 to —4.8 kms -1 in the blueshifted lobe and from 28.2 to 109.2 kms -1 in the 
redshifted lobe. The emission at lower velocities (|AV|<15 kms -1 with respect to the 
systemic velocity of 12 kms" 1 for MM2) suffers from side lobes and missing flux due to 
inadequate (u, v) coverage, hampering the interpretation of the morphology and kinematics, 
hence we focus on the high velocity CO emission in this work. The most remarkable feature 
in Figure [2] is a jet-like outflow emanating from MM2 to the southeast (SE), which exhibits 
extremely high velocity (EHV) emission up to —85.2 kms -1 for the blueshifted lobe and 
103.2 kms -1 for the redshifted lobe above a 3a level of ~60 mJybeam -1 per 1.2 kms" 1 . 
This outflow shows a knotty appearance with multiple clumps aligned along the flow axis. 
In addition, two more collimated flows, the northwest and northeast flows, are seen in Figure 
|2fe; the former apparently originates from MC, whereas it is difficult with the existing data 
to unambiguously assign a driving source to the latter. These two flows both extend well 
beyond the SMA primary beam of ~54" at 230 GHz, thus are not further discussed in this 
work. 

Figure [3] shows the SO (65-54) emission integrated over five velocity intervals. The 
blueshifted emission arising from the SE outflow is detected as an EHV clump in Figure [3^ 
and a low to high velocity clump in Figure [3b. In Figure [3b the blueshifted low velocity 
emission mostly traces a wide-angle outflow surrounding the radio jet axis. In Figure EU the 
redshifted low velocity emission reveals two compact clumps coincident with MM1 and MC, 
as well as extended structures probably attributed to outflows from MM1 and/or MC. In 
Figure [3k, a redshifted clump is detected at very high velocities; this structure can also be 
identified in CO in Figures [2k and[2J, and is apparently associated with MC. 



4. Discussion 

4.1. Physical and Chemical Diversities among the Dense Cores 

The CH3CN (12^-llft-) emission is a good thermometer of dense gas surrounding high- 
mass protostars. Assuming that all the K components are in local thermodynamical equi- 
librium and tracing the same gas, we perform a simultaneous fitting, taking into account the 
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optical depth effect, to all the detected components by means of grid search x 2 minimiza- 
tion. From Figure HI the best fitting models agree well with the observations, and yield gas 
temperatures of HOjJg K for MM1, 35j^j K for MM2, and 45±fg K for MC. 

In addition to the CH3OH (8_i,8-7o,7) E line shown in Figure [Tb, thermal emission in 
CH3OH (8o,8 — 7i ) e)J5 and (3_ 2 , 2 -4_i,4).E is detected in all the three cores but appears much 
brighter in MM1. The CH 3 OH (15^11-163^3)^ emission, however, is seen in MM1 alone 
(Figure [Tg). This is not surprising since this transition has an upper energy level (E up ) of 
374 K, dramatically higher than 39-96 K for the other three transitions. As indicated by 
the higher temperature inferred from the CH3CN emission, only the MM1 core possesses 
sufficient hot gas to populate the high-excitation CH3OH. 

HNCO shows a diversity similar to CH 3 OH; the HNCO (10 , 10^0,9) emission is seen 
in both MM1 and MC, and probably MM2 as well, but the HNCO (10i, 9 -9i, 8 ) emission is 
detected in MM1 alone. Lines of HNCO (JK a ,K b ) are suggested as a good probe of the far- 
infrared field since its K „.>0 ladders are popula t ed primarily by far-inf rared radiation rather 
than collisions with H 2 (IChurchwell et al.lll986l ; IZinchenko et al.ll2000l ). The warm/hot dust 
in MM1 is expected to produce far-infrared photons responsible for the population of the 
K a =l HNCO, giving rise to the HNCO (10i t g-9i t 8) emission, however, the cooler dust in 
MM2 and MC would not provide sufficient far-infrared radiation to populate appreciable 
HNCO at K a >0. 

The emission from SO2 and its isotopologues exhibits a significant diversity among 
the three cores. Besides the SO2 (ll5,7-12 4 , 8 ) emission shown in Figure QJ, the lines of 
S0 2 (22 7 ,i 5 -23 6 ,i 6 ), 34 S0 2 (111,11-100,10), 34 S0 2 (4 2 ,2-3i, 3 ), 34 S0 2 (22 2 , 20 -22i, 2 i), and 33 S0 2 
(lli^i-lOcio), are clearly detected in MM1. However, none of them is seen in MC or MM2. 
The detected transitions have E up ranging from 19 to 353 K, hence the temperature variation 
among the cores cannot solely account for the non-detection of any transition in MC or MM2. 
In addition, the SO2 chemistry does not seem to be much dependent on luminosity or mass 
of the central forming star, as bright SQ 2 emission was detected in low-mass star-forming 
regions as well (e.g., IRAS 16293-2422. iBlake et al.lll994l ). It is therefore the chemical effect 
that plays a central role in shaping the appeara nce of SOo in th ese c ores. If SO 2 form s 
via the reaction of SO+OH^S0 2 , as adopted by ICharnley I (119971 ) and iDoty et al.l (120021 ). 
the non-detection of SO2 points to the deficiency of OH in MC and MM2, since bright SO 
emission is detected toward all the three cores. This is plausible considering that the UV 



induc ed photodissociation of H 2 may be the primary source of OH in MM1 (ITappe et al. 



20081 ) and the formation of other molecules detected in MC and MM2 does not necessarily 
involve the OH radical. The rich emission lines from SO2 and its isotopologues seem to 
probe a more evolved stage of MM1, in agreement with recent sulphur chemistry models 
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f lDotv et alJl200a Beuther et al.1 120091 ) . 



With a gas temperature of 110 K and moderately rich molecular line emission, MM1 
seems to be at an early hot core phase (IKurtz et al.ll2000l ). Of the two younger cores MM2 and 
MC, the latter is likely more evolved with its more active molecular line emission. One issue 
to be addressed toward MC is that the core appears relatively faint in the dust continuum, 
which may be caused by the deficiency of dust or by a flattened density profile hence a 
large fraction of the emission being filtered out by the interferometer. With the existing 
data we cannot shed more light on this issue. However, its hot-core-like line emission is very 
likely excited by a high-mass protostar. Although nearby low-mass proto stars my excite 



emission in complex organ ic molecules such as CH3OH and CH3CN (e.g. iBottinelli et al. 



2004 ; IChandler et al.ll2005l ). the emission may not be detectable at much larger distances. 



Quantitatively, if the dust and gas reach thermal equili brium, the luminosity of the central 



sourc e can be estimated based on the heating of the dust (jScoville fc Kwanlll976l ; IZhang et al. 



20071 ). For dust temperatures of 35-110 K at scales of 1" and a typical dust emissivity index 
of 1.5, the luminosity of the central star ranges from 2xl0 3 to 4xl0 4 L G . Thus, the central 
protostars in the three cores are likely to be massive. 



4.2. The EHV Jet-like Outflow and Molecular Bullets 

To obtain a more reasonable assessment of the mass and energetics of the SE outflow, 
we calculate the gas mass down to |Ay|~6 kms -1 ; but to alleviate contamination from 
ambient gas and unrelated features we only take i nto account an ar ea closely encompassing 



the jet-like structure. Following the procedure of lOiu et al.l (120091 ). we derive the outflow 
mass, momentum, dynamical timescale, mass rate, and momentum rate of 0.22 M Q , 4.9 M Q 
kms -1 , 2.2xl0 3 yr, 10~ 4 M Q yr _1 , and 2.2xl0 -3 M kms -1 yr -1 , respectively. Considering 
the missing flux and optical depth effect at low velocities, the derived mass and energetics 
may represent lower limits. The overlapping blue- and redshifted emission along the jet-like 
structure suggests that the flow axis is strongly inclined to the plane of sky. In projection 
the flow axis appears to wiggle over a small angle of <20°. If the flow axis wiggles along the 
line of sight as it does in the plane of sky, the inclination angle is likely about 10°, half of 
the wiggle angle. The outflow mass and momentum rates corrected for the inclination angle 
would amount to 6xl0 -4 M Q yr _1 and 7xl0~ 2 M kms _1 yr _1 , respectively. The one-sided 
appearance of the outflow is probably ascribed to the peculiar environment. In projection 
the central source MM2 is close (~4") to the HH 80-81 radio jet, which is surrounded by 
a biconical cavity. The northwestern side of MM2 is likely deficient in dense gas and the 
outflow is expanding into the radio jet cavity, thus cannot create appreciable CO emission. 
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Figures[5K and[5b show, respectively, the mass- velocity (MV) and position-velocity (PV, 
cut along the major axis) diagrams of the SE outflow. The most interesting characteristics 
in the two diagrams are the distinct EHV features that appear as abrupt bumps in the MV 
diagram (denoted by arrows) and as enhanced condensations in the PV diagram (denoted by 
dashed circles). This kind of compact and fast mo ying features are recognized as "molecular 



bullets" in low-m ass class protostellar outflows (IBachiller et al.lll990l ; lHatchell et al.lll999 



Nisini et al.l 120071 ). The bullets in the SE outflow are, to our knowledge, the first discovery 
of molecular bullets originating from a high-mass star-forming core. The bright Bl and 
B2 bullets have sizes of ~0. 04-0. 05 pc and masses of 10~ 3 M . However, recent multi- 
transition SiO and CO observations yi eld hot gas (>100 K) in low-mass molecular bullets 
(lHatchell et al.lll999l ; iNisini et al.l 120071 ) : if this is the case for the Bl and B2 bullets, their 
masses would reach >10~ 2 M . It is also worth noting that the Bl bullet is clearly detected 
in SO (65-54), suggesting shock activity in this bullet (Figure Eh)- 

From its velocity profile shown in Figure [5b, the B2 bullet has the bulk material moving 
at highest velocities, with a gradual decrease toward the systemic velocity. The Bl bullet 
exhibits a similar but less prominent profile. The velocity structure of the two bullets is 
in direct contrast to what is expected if the bullets were created in situ from entrained or 



Zhane et al. 


1995: 
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2007) 



the systemic velocity (e.g. 
majority of the bullet material is ejected from the close vicinity of the central protostar. In 
this scenario the presence of multiple molecular bullets directly manifests an episodic nature 
of the outflow. 



Marti et al.l (119991 ) resolved H2O maser spots toward MM2 into two groups, with one 
group aligned in a linear distribution and the other located to the NE. Although the H 2 
maser spots are seen very close to MM2 (~500 AU in projection), we find that in both veloc- 
ity and orientation the NE and linear maser groups agree well with the Bl and B2 bullets, 
respectively. This further supports the scenario that the bullets are ejected from the close 
vicinity of the central protostar. A sketchy and tempting interpretation is that the episodic 
ejections from the central protostar or star-disk system collides with the existing flow, re- 
accelerating the flow material to form the fast moving bullets; the maser emission, which is 
pumped under critical conditions, signifies the collision sites. Future proper motion mea- 
surements of the H2O maser spots and multi-transition observations of the bullets may place 
more stringent constraints on the excitation mechanisms of the two extreme phenomena. 

High velocity, jet-like molecular outflows are often seen in low-mass class protostar s 



{ e.g. , iGuilloteau et al.lll992l ; IBachiller et alJll995t iGueth fc Guilloteaulll999l : iLee et al.ll2007[ ) ; 
in theoretical models they are physically linked to the earliest-phase disk accretion within a 
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dense envelope ( jShang et al.ll2006l ; lFendtll2009l ). However, jet-like outflows with hi gh velocity 



emiss ion are still rarely identified toward high - mass star-forming regions (e. g., iQiu et al. 



20071 ). Based on single-dish CO observations, iRodriguez-Franco et al.l (119991 ) report high 



velocity bullets and jets toward the well-known Orion-KL outflow, but rec ent interferometric 
observations show a fan -shaped or explosive morphology of the outflow (IBeuther fe Nissen 



20081 ; IZapata et al.ll2009l ). The SE outflow has a well-defined jet-like morphology, reaching an 



order of 100 km s _1 in both blue- and redshifted lobes, thus apparently favors a disk-mediated 
accretion onto the central protostar. The EHV emission may trace part of the outer layer of 
an underlying wind. Moreover, the outflow contains multiple molecular bullets, indicating an 
episodic mass-loss process and further pointing to an episodic nature of the central accretion 
process. 



5. Summary 

We report the discovery of jet-like CO outflows in the HH 80-81 star- forming region. In 
particular, the outflow arising from core MM2 exhibits blue- and redshifted line wings close to 
100 kms -1 from the cloud velocity, and molecular bullets, which points to an episodic, disk- 
mediated accretion for massive star formation. Three dense cores identified in continuum 
and molecular line emission demonstrate physical and chemical diversities which suggest an 
evolutionary sequence that MM2 is younger than MC, which is in turn younger than the hot 
molecular core MM1. 

We thank Ciriaco Goddi for discussions on water and methanol maser emission in high- 
mass star-forming regions. 
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Fig. 1. — (a) The 1.3 mm continuum emission, contouring at (—1, 1, 2, 4, 7, 11, 16, 22, 29, 
37, 46, 56, 67, 79) x 0.006 Jybeam -1 ; a filled square denotes the brightest H 2 maser spot 
from Marti et al. (1999); a plus sign marks the position of the central radio source from 
Marti et al. (1999). (b)-(i) The molecular line emission, integrated from 9.6 to 25.2 kms -1 
and contouring at (-1, 1, 2, 4, 7, 11, 16, 22, 29, 37, 46, 56, 67)x0.26 Jybeam" 1 kms" 1 for 
CH 3 OH (8-1^-70,7) E, and integrated from 10.8 to 16.8 kms -1 and contouring at (—1, 1, 2, 3, 
...)x0.12 Jybeam -1 kms - 1 for the othe r lines ; a filled triangle in panel (b) denotes a 44 GHz 
methanol maser spot from lKurtz et al.l (120041 ); the line transition is labelled in the upper left 
of each panel. Hereafter the three stars denote the peak positions of the three cores and a 
filled ellipse in the lower left delineates the FWHM of the corresponding synthesized beam. 
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Fig. 2. — The CO (2-1) emission integrated over every 20.4 kms -1 , as indicated in the 
upper left of each panel, from —85.2 to —4.8 kms -1 and from 28.8 to 109.2 kms -1 ; the first 
and spacing contours are 0.4, 0.6, 1.2, 1.8 Jy beam -1 kms -1 for panels (a) and (h), (b) and 
(g), (c) and (f), (d) and (e), respectively; hereafter the double dashed arrows delineate the 
orientation of a thermal radio jet from Marti et al. (1993); the large dashed circles mark the 
primary beam of the SMA. 
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Fig. 3. — The SO (65-54) emission, integrated over five velocity intervals as indicated in 
the lower left of each panel; the first and spacing contours are 0.24, 0.3, 1.2, 0.3, 0.4 
Jy beam -1 kms -1 for panels (a) to (e), respectively. 
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Fig. 4.— The CH 3 CN (12-11) spectra for MC (top), MM2 (middle), and MM1 (bottom), 
taken from the positions indicated by the three stars in Figure [U The best fitting models 
are overlaid in thick gray lines. 
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Fig. 5. — (a) Mass- velocity diagram of the SE CO outflow, calculated in each 1.2 kms -1 
channel; gray dots and dark squares denote measurements of the blue- and redshifted lobes, 
respectively; three arrows mark enhancements attributed to the molecular bullets, (b) 
Position-velocity diagram along the major axis; three dashed circles highlight the molec- 
ular bullets; the vertical dashed line denotes the systemic velocity of the MM2 core, (c) The 
velocity profile of the B2 bullet. 



